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Threshold current of domain wall motion under spin-polarized electric current in ferromag- 
nets is theoretically studied based on the equation of motion of a wall in terms of collective 
^ ' coordinates. Effects of non-adiabaticity and a so-called /?-term in Landau-Lifshitz equation, 

which are described by the same term in the equation of motion of a wall, are taken into 

f^ ' characterized by different dependence of threshold on extrinsic pinning, hard-axis magnetic 

o 



account as well as extrinsic pinning. It is demonstrated that there are four different regimes 



anisotropy, non-adiabaticity and /3. 
j^ . KEYWORDS: spin torque, domain wall, threshold current, MRAM 



1. Introduction 

Intensive studies, both theoretical and experimental, have been carried out to understand 
and control current-driven domain wall dynamics in nanoscale magnets. ^"^^ So far all the 
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^ experimental results ' ' indicate that wall motion can really be induced by a current density 

larger than the threshold, jci and this threshold is of the order of 10^^[A/m^] in metals^' ^"^ 
and of 10^[A/m^] in magnetic semiconductor.^ The key issue yet to be understood is the origin 
of this threshold and driving mechanism. 

Theoretically, current-driven domain wall motion was predicted and studied by Berger 
since 1978.^"^^ He showed^^ that the motion in the limit of thick wall (adiabatic limit) is 
induced by the transfer of spin angular momentum from conduction electron to the wall via 
exchange interaction (spin transfer). Another driving mechanism due to electron scattering 
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(momentum transfer) was also pointed out on a phenomenological ground.^'' Wall dynamics 
was discussed in terms of two variables, wall position X and polarization (angle from the easy 
plane) (j), namely, within a rigid wall approximation introduced by Slonczewski.^^'^^ Recently, 
this problem was reformulated from a microscopic point of view by two of the authors^^ 
within a rigid-wall approximation. The work of Berger was extended there to systematically 
incorporate non-adiabaticity, and the momentum transfer effect was shown to be governed by 
the wall resistance, R^. A solution in the adiabatic limit in the absence of extrinsic pinning 
was presented there and existence of threshold current was pointed out.^^ This threshold is 
due to an intrinsic pinning arising from a deformation (change of 0), and the threshold was 
shown to be given by hard axis anisotropy (or demagnetization field), K^, as 

^ = ^"■-^' « 

where A is wall thickness, P = js/j is the polarization of current, S and a being the magnitude 
of local spin and lattice constant, respectively. 

Recently a new torque term in Landau-Lifshitz-Gilbert (LLG) equation was proposed by 
Zhang and Li"*^^ and Thiaville et al.^'^ This torque, proportional to [3S x (jg • V)S, where jg is 
spin current and /3 is a coefficient, is perpendicular to the spin transfer torque, (jg • V)S', and 
is argued to arise from spin relaxation of conduction electron, ^^'^^ modification of damping by 
current, ^^ etc.^*^ This term (sometimes called a (3 term) was shown-*^^'^^ to remove the thresh- 
old of the intrinsic origin, and that the wall starts to move for any small value of spin current 
with a finite velocity given by X oc — jg, where a is Gilbert damping parameter, if extrinsic 
pinning is absent. Smearing of intrinsic threshold occurs also due to non-adiabaticity. This is 
seen from the fact that (3 term plays exactly the same role as momentum-transfer effect on the 
domain wall dynamics described by X and (p (i.e., as far as wall deformation besides (/>-mode is 
neglected). ^'^ Existence of such term was shown theoretically in a different system of ferromag- 
netic j unction. ^"^ At finite temperature, we have to take account of thermally assisted process, 
which also contributes to smoothen the velocity-current curve near the intrinsic threshold, 
j],}'^ In refs.,^^'^^ value of (5 was discussed to be identical to a (in the case of a single band 
model in ref.^^). In the present paper, we neverthelss treat j3 and a as independent, since in 
reality, multiband effect^^ and contribution to a from non-electron origin would result in a 
difference. In the case of rigid domain wall we consider, [3 is effectively replaced by (3' below, 
which includes the effect of momentum transfer, and thus (3' and a should be regarded as 
independent. 

Once non-adiabaticity or the (3 term are taken into account, pinning of extrinsic origin 
(such as defects) becomes essential in the dynamics, since extrinsic pinning blocks the motion 
at low current. This was indicated in a simulation by Thiaville et al.^ They simulated extrinsic 
pinning by surface roughness of the wire and showed that threshold current depends on 
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roughness. Ohe and Kramer^^ recently studied a wall motion solving fully the torque arising 
from the conduction electron numerically. Non-adiabaticity (momentum transfer force) was 
thus fully taken into account there. Extrinsic pinning being neglected, the threshold obtained 
was zero, consistent with the result of Thiaville et al, where the non-adiabaticity was effectively 
taken account of in terms of /?.^^ These analysis so far are not enough to draw conclusion on 
the relation among threshold, /3-term and extrinsic pinning. The aim of the present paper is 
to explore the wall dynamics with both (3 term (and non-adiabaticity) and extrinsic pinning 
taken into account on a basis of the equation of motion of a domain wall. 

In this paper, the sum of /3 contribution and non-adiabaticity, i.e., the coefficient /3' (Eq. 
(4)), is assumed as positive. Positive /?' seems reasonable at present in that all the so far 
proposed origins of /3, spin relaxation^^ and modification of damping, ^^'^^ leads to positive P 
and the effect of non-adiabaticity (below Eq. (4)) is also positive. Some of the results in the 
weakly pinned regime are modified if /3' < as we mention below. 

2. Equation of motion of domain wall 

The equation of motion of local spin under current is written as^^' 

S = B,sXS+^SxS-^{js-V)S-^[Sx (js • V)S] + T„a + Tpin, (2) 

where B^s is the effective field arising from spin Hamiltonian, a represents Gilbert damping, 
which is introduced phenomenologically. Spin torque (spin-transfer effect) is represented by 
a term (js • V)S^. This term is derived in the adiabatic limit, ^^ and thus the spin S in the 
above equation of motion must be slowly varying compared with Fermi wavelength, kp . Tna. 
and Tpin represent torque from the non-adiabaticity and pinning, respectively. Explicit forms 
of these two terms are discussed when the equation of motion of domain wall is discussed. In 
fact, non-adiabaticity is not easily expressed in the equation of motion of a spin, since Tna is 
non-local in space, "^^ while its effect on collective object such as domain wall is clearer as was 
discussed in ref.^^ 

We consider a planar domain wall and take account of only a deformation mode described 
by the polarization angle (f> out of the easy plane. This is justified if the wire is narrow (L± < A, 
L± being the wire width) and if the hard-axis anisotropy, K±, is weaker than the easy-axis 
one, K.^'^ Such wall is described by the collective coordinates, X and 0, and the equation of 
motion derived from Eq. (2) reads^^' ^^' ^'^ 

Q 

Pj 

(3) 

where Vc = K±XS/{2h) corresponds to the drift velocity of the electron spin at j = jl, and 
P = js/j is the polarization of the current. Left-hand side of the second equation represents 
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total force acting on the wall. The force contribution from the current j is a sum of /3- 
contribution and non-adiabatic contribution (in other words, momentum-transfer effect);^^ 

P' = /3P + /3na, (4) 

where Pna = 2?i5^^"''^"^ ^^ ^ dimensionless measure of wall resistivity, n is the electron 
density, / = jA, A being the crosssectional area of the wire. Pinning potential is treated as 
harmonic, with frequency of 0, and range S,', 

V^pin = ^M^n^X^ - eWi^ - \X\), (5) 

where M„ = ^ , 3 = ^ -^ is wall mass, N = 2XA/a'^ is the number of spins in the wall. 
In most cases, pinning range ^ is of order of width, ^ ~ A. Pinning term in Eq. (3), /pin, is 
defined as /pin = -pj^-Fpin with Fpin = — qx^^ ^^'^ given as 

The wall velocity in the limit of large current is obtained from Eq. (3) as 

X I a^ f P 



7T7^ + a/3' i 



A 1 + q2 eA V 25 

; I a? f ^, P 



/3' - «7T7^ J- (7) 



l + a2eAV 25 

Thus the spin torque efficiency, defined as i] = X -fp. , can be larger than unity when P' 
contribution is taken account. 

3. Numerical results 

In this section, results of numerical calculation are presented. The equation of motion 
solved is in terms of dimensionless parameters; 

— (^X-a<p) = sm2(l) + Pj 

^,(</. + ax) = -ln^xei^-\X\)+(3'j = -VoX9{^-\X\) + (3'j, (8) 



where t = tvJX, X = X/X, Q = QX/v^ = ^f V^' ^ = ^' ■^' = ^J ^^^ ^0 = ^^ 
Vo^^ = ^^ {l)\ with Vo ^ ^n-e = ^¥z {if = f €? being a pinning poten- 
tial strength per spin. Calculation is done with a = 0.01, S = 0.5^^ and P = 1. Threshold 
current is plotted as function of 17 and Vq in Figs. 1 and 2 for several values of (3' . It is clearly 
seen in Fig. 1(a) that behaviors for /?' = and /?' 7^ are quite different except for extremely 
strong pinning regime (Vq > 100 ~ 1/a). We will come back to this difference later. 
Let us focus on the case /?' 7^ first. Fig. 1 indicates that there are three regimes; 
I) Weak pinning regime ((7 < 0(1)), ic grows with Vq 
II) Intermediate regime (0(1) ^ ^ ^ 0(a^^'^), jc ~constant ~ 0.7 — 0.8 
III) Strong pinning regime {Cl > 0{a~'^''^)), jc oc Vq 
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Fig. 1. (Color online) (a): Threshold current jc plotted as function of pinning frequency Jl = •\/2yo 
for several values of /?' with a — 0.01 and P = 1. (b): Threshold current in weak pinning regime. 
Fitted curves for /3' > 0.1 are jc oc Cl^. For small /3'(< 0.02), jc is linear in Cl. 





Fig. 2. (Color online) (a): Threshold current jc in weak pinning regime plotted as function of pinning 

Vg 
2(3' 



strength Vq for several values of /?' with a — 0.01 and P = 1. For /3' > 0.1, jc — -^ is linear in 



Vq. (b): Weaker pinning regime. For small /3'(< 0.02), jc oc yVo. 

It should be noted here in figures that (3' affects the threshold only for weak pinning, Vq ^ 1 
(i.e., Vq < K±). Looking closer, we see that the weak pinning regime actually consists of two 
regimes; 

I-a) Small /3' /Q: jc oc ^A^ 

I-b) Large (3'/^: ^ oc Vo/(3' 

Physical picture of these regimes is as follows. In regime I, current is low, so </> does not 
grow, and dynamics is described by X. In regime I-a), /?' is negligible and depinning is due to 
spin-transfer. Namely, spin-transfer term (proportional to P) gives a finite velocity to domain 
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wall and the wall escapes from the pinning potential by use of initial kinetic energy supplied 
from the spin-transfer torque. In other words depinning is due to a kinetic energy supplied 
by spin transfer. In regime I-b), depinning is governed by a force. It occurs when the force 
due to current exceeds pinning force, as is seen from the expression jc oc Vo//3'. Dynamics in 
regime II) and III) is described by (p. Regime II) is a regime where depinning is due to spin 
torque, but terminal velocity is determined by /?'. Depinning mechanism here is the same as 
the intrinsic pinning pointed out in ref.,^^ but the lowering of jc by a factor of 0.7 ~ 0.8 by 
P'-terui is one of new results of the present paper. This is the reason why jc looks different 
for /3' = and /3' 7^ in Fig. 1 (a). Regime III) is a strong pinning regime argued in ref.^^ 
Depinning here is due to spin-transfer torque. In both regimes II) and III), /?' does not affect 
the threshold current. 

Let us look into each regime closely. 

4. Analytical results 

4-1 Weak pinning regime : I) 

Under weak current, J ^ 1, </> remains small and wall dynamics is well described by X 
only. This is the first regime I). Linearizing sine-term in Eq. (3) as sin2(^ ~ 20, can be 
eliminated to obtain a simple equation for X as^'^^ 

{l + a^)d'lX + \diX + ^^X = F-^, (9) 

where 1/f = 2a ( 1 + ^Vq J and 

Fj ^ 201 (10) 

is a dimensionless force due to current. We consider a case of steady current and weak damping; 
2rif > 1 A general solution to the equation (9) (inside the pinning potential) is given as 

X{t) = (U- + e'"^ {Acosn'i + B sinn'i), (11) 

$7^ 



where O' = \/^'^ — j^ and A, B are constants. Initial condition required is A(0) = and 
diX{Q) = Pj. The second condition on the wall speed comes from the first equation in Eq. 
(8) (with (p{0) = 0), and is the most important consequence of spin-transfer torque; namely, 
spin-transfer torque gives initial speed to the wall. With these initial conditions, we obtain 

X(t) = ^ (l-e'2T (cosn'i + -L^ smn'i)] + Se"^ sinll't. (12) 

The first part is governed by a force from /?' and the second is driven by a spin-transfer torque 
term. In the case of small [3' , spin torque contribution leads to a maxmum displacement 

l^maxl - £j, (13) 
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while 

l^maxl ^ -^j, (14) 

if /?' is large. (We assumed here that damping is weak (fiV ^ 1).) The first regime corresponds 
to regime I-a) and the second to I-b). Threshold current in each case is given as 



•la) ^ ? 

■Jc py 

and 



(15) 



4/3' A 
The crossover occurs at 



jI'' - Si (16) 



2 
In terms of dimensionful quantities, 



P'c ^ ?^- (17) 



I,) ^ 2V2S e VKIV, _ 2V2JW, .^g. 

and 

. lb) ^ ^^^. ^ SeVo_^^ ^ P±V^X, 
■'" 4.v^a?Wr ha^ \P'\ i 5|/3'|Kxe 

Note that simple comparison of pinning force and F-^ in Eq. (9) gives a result correct up to a 
numerical factor, jj^^ ^ i^^^. 

The pinning strength Vq is experimentally accessible by driving the wall by magnetic field. 
Magnetic field B along easy axis add a term in Eq. (3) 

fs = '-^B, (20) 

where g = 2. By a simple comparison of pinning force and magnetic field, Vq is written in 
terms of the depinning magnetic field Be as 

Vo = ^gf^BBcj, (21) 



and so jj^' is simplified to be 



^"" = I^Y'^-^^'W (^^' 



4-2 Intermediate regime : II (Intrinsic pinning) 

This regime could be important for application since the threshold is not sensitive to the 
sample irregularities. Depinning in this regime j > 0(1) is described by <j) as done in ref.^^ 
The reason is that the effective mass of (/>-" particle" , given by l/Vo^^(see Eq. (23)), becomes 
lighter than the corresponding mass of X-" particle" given by l/i^_L, and so 0-" particle" is a 
better variable to describe dynamics for strong pinning. By eliminating X from Eqs. (8), we 
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obtain 

(1 + a2)a?0 + ad^tj) (2 cos 2<j) + Vq) + Vq sin 2(j) + jVqP = 0. (23) 

Thus /?' does not affect the dynamics of (p. (Correctly, this feature is specific to a harmonic 
pinning potential, and anharmonicity results in appearance of /3'. In fact, /3'-term is eliminated 
from the equation of motion if one replaces X in Eq. (8) by X' = X — -£^j (i.e., shift of stable 
point of X). Even in anharmonic case, nevertheless, we have numerically checked that the f3' 
does not lead to important contribution in this regime.) 

Prom Eq. (23), we see that the energy barrier for (p vanishes when jc ~ P~^, irrespective 
of pinning strength. Once (j) escapes from local minimum, its velocity is given by Eq. (23) as 

di4> ^ — • (24) 

a 

This corresponds by use of Eq. (8) to a maxmum displacement of the wall of 

Xm..^-^di<P^^. (25) 

For intermediate pinning strength, aVo < 1, iXmaxI exceeds ^/X (if ^/A ~ 1), i.e., depinning of 
X occurs as soon as (p is depinned. Thus the threshold is roughly given by jc ~ P ^, but this 
estimate turns out to be too rough to determine numerical factor correctly. In fact, numerical 
result in Fig. 1 indicates that jc is actually given by 

ie~0.7xP-i (;3'^o), (26) 

while jc ^ P^^ when (3' = 0. Result of /3' = case is in agreement with ref."*^^ The reason 
for the difference between /3' = and /?' 7^ is understood as follows. In the case of /?' = 0, 
even if X escapes from the pinning center at current j > 0.7 P^^, terminal velocity vanishes 
if j < P~^, since the motion stops due to the intrinsic pinning effect (i.e., (j) reaches a steady 
value and X becomes zero). On the other hand, if /3' 7^ 0, steady motion of X is possible as 
soon as X escapes from the pinning (j > 0.7P~^). This is the reason the threshold value in 
intermediate regime is different for /?' = and /3' 7^ (Fig. 1). 

If P' < 0, or more precisely, if the relative sign between /?' and P in Eq.(8) is negative, 
the /3'-term will drive the depinned wall back to the pinning center, and the threshold in this 
regime is given by the intreinsic value jc = 1. 

4-3 Strong pinning regime : III 

Eq. (25) indicates that for extremely strong pinning, Vq > a^^, the wall is not always 
depinned even after (j) escapes from the potential minimum. Depinning occurs at 

Jc ~ -^^, (27) 



as has been pointed out in ref. 
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5. Wall speed 

5.1 Close to threshold 

Another important quantity to understand mechanism of waU motion is the wall velocity 
at jc- We consider a steady current or a pulse with duration long enough (longer than time 
scale of r), and the weak pinning case, jc ^ j\. (i.e., regime I). The terminal wall velocity 
after depinning (j > jc) is obtained from eq. (9) with fi = as 

V ^ diX = ^j. (28) 

id! 

At jc, the wall velocity suddenly jumps from zero to — jc, and so the jump in dimensionless 
unit is given by Ati^^-* = ^-^^ in regime I-a) and A-u^^' = — — ^ in regime I-b). In terms of 
dimensionful quantity, Av/jc = ^^^ ^^ have 



Av'-) = ^' ^^SVKlVo ^ 



a P h ' 

^^ib) ^ i^^. (29) 

ant. 

These behaviors are seen in Fig. 3. Comparing with result of ref.,^''' behavior appears qualita- 
tively similar in the large (3' regime. Detailed comparison is, however, not possible since value 
of K_\_ and Vq in ref.^^ are not known. It should be noted that the wall velocity near extrinsic 
threshold is discontinuous at zero temperature. This is because the wall as soon as depinned 
feels a tilted potential (due to (3') and has finite velocity. This might explain rather scattered 
experimental data for velocities in metals.^' ^ In contrast, results on semiconductor indicates 
quite smooth temperature dependence,® which might suggest different origins of threshold 
from those in metals. 

5.2 After depinning (or no pinning) 

Let us see how wall velocity behaves after the wall is depinned. We solve here Eq. (8) with 
Vo = 0, which is a simple equation of given as 

d-^cj) = ({(3' - Pa) J - a sin 2(p) . (30) 

For a steady current, this equation has a solution of 



sm2(P= {—-P]j-sgn(P-—j) 5L±^^ ^J (31) 

*a yj B V a'^\{f3' -aP)j\+asmi2u;t-^) ^ ' 



p-i\f 



where uo = , " ■> 

velocity is given from Eq. (8) as 

1 

1 + q2 



1, sint? = ^ - "^^, -^ , and sgn(i?) denotes sign of B. The wall 



|P-^il • 



d^X = rr—p^ (sin 20 + (P + a^')j) . (32) 



9/15 



J. Phys. Soc. Jpn. 



Full Paper 



We see that anomaly in the velocity appears when uj switches from real to imaginary, i.e., at 

j = ja, where 

1 



Ja 



P 



/3' 



(33) 



Above and below jg,, wall dynamics is quite different. Above jg,, wall velocity (32) has an 
oscillating component, while wall reaches a steady motion if below j^. The average velocity 
for j > Ja is obtained by use of Eqs. (31) (32) and 

1 r . 1 1 



rT 
— I dt 
'o 



T./n ^Asm(2Tr^)+B ^B'^ - A^ ' 



T) 



(for \B\ > ^ and T = 2t:/uj) to be 



a 1 + Q^ 



P' 



)J] 



P 



H' 



a 



1. 



sin2(^ 



PVi. 



Below jai sin2(/) approaches at t ^ cx) a steady value of 

and so terminal velocity is obtained as 

a 



(34) 



(35) 



(36) 



(37) 



Note that the analysis here is the case of Vq = 0. In the presence of extrinsic pinning, velocity 
vanishes below jc as seen in Fig. 3. 



The anomaly at ja ^-re seen in Fig. 3 in the case of /3' = 0.001 (at j 



0.9 



1.1) and 



/?' = 0.02 (at i = 1), and the behavior is in agreement with numerical result of ref.,"*^^ where 
the anomaly of velocity under current was first reported. This anomaly is essentially the same 
as the anomaly under magnetic field, known as Walker breakdown. ^^'^^ 

Let us here briefly consider the effect of external magnetic field along easy axis. As seen 
from in Eq. (20), the field replaces I3'~j by /3'i + 6, where h = -^^^fjf^ = -^f^ (^ = 2). 
The average wall velocity with Vq = is then given if j > ja by 



{X) 



P'j + b , sgn[[(P-^)j]-A] 



+ 



a 



and by 



l + a2 



{X) 



P 



H' 



a 



1, 



/3'j + & 



a 



if j < Ja, where the anomaly now occurs at 



the lb denotes sgn[(PQ — /3') — 6/ja]- 



b± a 
Pa-P'' 



(38) 



(39) 



(40) 
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Fig. 3. (Color online) Wall velocity as function of current for il — 0.5 and a — 0.01. A jump in 
wall velocity is seen at j — jc- Crossover from l-a){jc ~ ^//3') to I-b)(jc ~ (l) regime is seen at 
/3'~ §-0.1. 



Table I. Summary of threshold current. The second to the last column indicates the mechanism of 
depinning, either spin transfer or force {/3'-term). The last column is a "good" variable to describe 
depinning. We see that purely intrinsic pinning occurs in regime II, while other regimes are affected 
by extrinsic pinning. 



regime 




threshold 


depinning mechanism 




I-a: Weak pinning 


n<o{i),(3'<o{n) 




spin transfer 


X 


jc oc \/-f^±Vb 


I-b: Weak pinning 


n<o{i),p'>o{n) 


Jc oc Vo/P' 


/?' 


X 


II: Intermediate pinning 


o{i)<n<o{a-^) 


jc oc K± 


spin transfer 


^ 


III: Strong pinning 


0{a-')<Q 


jc oc Vo/a 


spin transfer 


<P 



6. Discussion 

Let us summarize the results in table I. It is interesting that such a simple set of equation 
of motion results in so rich behaviors. This fact is very important for device application, since 
reduction of threshold current, which is a must to develop MRAM based on current-driven 
domain wall motion, can be realized in different ways depending on system. For instance, in 
regime II), threshold is governed by the average sample shape, and not by the sample quality 
(roughness) as pointed out in ref.^^ In contrast, region I-b) is most strongly affected since 
magnetic imperfection results in Vq, and magnetic defects modifies P' strongly by affecting 
the electron transport. Regime I-a) is moderately affected by both sample shape {K±) and 
defects (Vq). Existence of these different threshold may explain the controversy among exper- 
imental results. In fact, some experiments indicate threshold is very insensitive to artificially 
introduced pinning centers, ^^ while some experiments indicates pinning is essential. ^^ 
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So far experimental results on metallic samples all showed that threshold current is of order 
of 10"'^^[A/m^]. If we use experimentally estimated value of K_\_ ~ 0(1)[K],^^ the observed 
threshold is orders of magnitude (10^^ — 10^^ times) smaller than the intrinsic threshold, 
j],. For instance, a sample of Yamaguchi^' ^^ showed jc = 1 x 10"'^^[A/m^]. The anisotropy 
energy is estimated to be K^ = 2.4[K], and using S" ~ 2) ^ — 2.2A and P ~ 0(1)) we 
obtain j\. = 5.8 x 10"'^^[A/m^], i.e., jc/jc ~ 0.02. The observed low threshold in metals thus 
should be regarded as due to an extrinsic pinning in regime I-a) or I-b). Let us first try 
to explain experimental result^^ assuming regime I-a). We assume .^ ~ A. Pinning potential 
is estimated from the measured depinning field of B^ = 0.01 — 0.1[T]. By use of Eq. (21), 
Vq = 0.34 X (10-2 ~ 10-i)[K] = 4.7 x (lO^^^ _ l0-25)[j]^ i.e., ^ = 1.4 x (lO^^ _ iQ-2y 

Thus jc = (0.21 ~ 0.67) x j],. This value is still too big to explain the experimental value. 
Velocity jump becomes /S.v ' = ^ x839[m/s], so extremely small (3' (^ ~ 4x 10^^) is required 
to explain experimental value of At; ~ 3[m/s].^ If we assume regime I-b), the threshold is 
Jc^^-* = -njTT X 2.8 X (10^'^ ~ 10^2) xj].. Experimental value could be reproduced if /3' = 0.1 ~ 1. 
But such large value of (3' cannot be explained within the current understanding that fi' arises 
from either non-adiabaticity^^' -'^'^ or spin relaxation. ^^ Instead, Af cannot be explained by 
use of the above Vq assuming I-b), as it predicts too large value of /S.v ' = 10^[m/s]. Thus, 
honestly, none of the above predictions are successful in explaining experimental result of 
metals quantitatively. 

There remain some possibilities to resolve this disagreement. For instance, estimate of 
Vo by use of experimental Be could be an over-estimation if effective barrier height Vq is 
greatly reduced by heating under current, while such heating does not occur under static 
magnetic field. Let us estimate the pinning potential which gives the experimental value of 
jc- Assuming regime I-a), experimental value of jc/Jc ~ O-O^ is reproduced if ^ = -^ = 
1.3 X 10"^, which corresponds to Vb = 3 x 10-^[K]= 4.5 x 10~^[T]. This is two orders of 
magnitude smaller than the value extracted from B^. For I-b), we have fi = (3' x 10"^. From 
the experiment, Av/{^jc) = 3[m/s]/67[m/s]= 0.05. This value is equal, for regime I, to ^, 
so /?' = 5 X 10^'^ if a = 0.01. So in case I-b), fi = 5 x 10^®. Thus, assuming either regime I-a) 
or I-b), the experimental results could be explained by an extremely weak pinning potential, 

p- = 10-6 ^ 10-5. 

The situation is very different for magnetic semiconductor.^ The observed jc is in quan- 
titative agreement with intrinsic threshold, j'^, and the velocity above threshold is also in 
agreement with spin-transfer mechanism. Thus extrinsic pinning and /3' does not seem to play 
roles. 

Identification of the origin of threshold is thus necessary to lower the threshold current in 
experiments. It would be very important first to see if thresholds observed in metallic systems 
are more or less governed by Vq and (3' , and second to try to lower threshold by controlling 
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sample quality (Vq) and (3' by introducing heavy-atom impurities. 

The authors are grateful to Y. Yamaguchi, T. Ono, M. Yamanouchi, H. Ohno, M. Klaui, 
A. Thiaville, J. leda, N. Nagaosa, J. Inoue, S. Maekawa and S. Parkin for valuable discussion. 
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